New Class of Magnetoresistance Oscillations: Interaction of a Two-Dimensional 
Electron Gas with Leaky Interface Phonons 
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We report on a new class of magnetoresistance oscillations observed in a high-mobility two- 
dimensional electron gas (2DEG) in GaAs-AlrGai-^ As heterostructures. Appearing in a weak 
magnetic field (B < 0.3 T) and only in a narrow temperature range (2 K < T < 9 K), these 
oscillations are periodic in 1/B with a frequency proportional to the electron Fermi wave vector, 
kF- We interpret the effect as a magnetophonon resonance of the 2DEG with leaky interface-acoustic 
phonon modes carrying a wave vector q = 2kp- Calculations show a few branches of such modes on 
the GaAs-AUGai-^As interface, and their velocities are in quantitative agreement with the data. 

PACS numbers: 73.50.Rb, 68.35.Ja, 73.50.-h 



There are several classes of transverse magnetore- 
sistance (MR) oscillations known to exist in a two- 
dimensional homogeneous electron gas (2DEG). The 
most common of these are the Shubnikov-de Haas oscilla- 
tions (SdH), which arise from a magnetic field B-induced 
modulation of the density of states at the Fermi level Ep. 
They become more pronounced with decreasing temper- 
ature T. The magnetophonon resonance (MPR) ||, ^ is 
a source of another class of oscillations resulting from the 
absorption of bulk longitudinal optical phonons. These 
resonances appear under the condition lolo = l<-<Jc, where 
ujlo and lo c = eB/mc are the optical phonon and cy- 
clotron frequencies respectively, I is an integer, and m is 
the effective mass of the carriers. These oscillations are 
only seen at relatively high T ~ 100 — 180K gj. Both 
SdH and MPR are periodic in 1/B, but the SdH fre- 
quency (reciprocal period) scales with electron density 
as n e , whereas MPR is n e -independent. 

In this Letter, we report on a new class of MR oscil- 
lations H observed in a high-mobility 2DEG in GaAs- 
AlzGai-zAs heterostructures. Unrelated to either of the 
above origins, these novel oscillations are still periodic 
in 1/B, but they appear only in a narrow temperature 
range (2 K < T < 9 K), and their frequency scales with 
y/n^. We interpret the data in terms of a magnetophonon 
resonance mediated by thermally excited leaky interface- 
acoustic phonon (LIP) modes. In principle, the surface 
modes might provide a good explanation as well, but in 
our case 2DEG is located so far from the surface (~ 0.5 
yitm) that no such interaction is possible. 

The leaky interface modes have been studied for a few 
decades in connection with the Earth's crust The 
term "leaky" shows that the waves propagate at a small 
angle with the interface, so that the energy radiates away 
from the boundary. For some specific parameters these 
waves may not be leaky la], but for the interface under 



study all of them are leaky. Despite the fact that LIP 
is commonly present in layered material systems (6|, it 
has so far not been considered on the GaAs-ALjGai-^As 
interface. Due to radiation of energy, the frequency and 
velocity of leaky waves are complex: u — u> / q = ur — iuj 
with ui <C itij- 

The novel oscillations can be explained by a simple mo- 
mentum selection rule which is derived later in the paper. 
It states that at high Landau levels (LLs) the electrons 
interact predominantly with the interface phonons carry- 
ing a wave vector q — 2k p, where kp is the Fermi wave 
vector of the 2DEG at zero B field. The condition for 
resonant absorption or emission of an interface phonon 
is then given by 



2kpUR — Ilo c , I = 1,2, 3, 



(1) 



We claim that Eq. (|lj) determines the values of B for 
the maxima in these new MR oscillations. It shows that 
the oscillations are periodic in 1/B with a frequency / = 
2kFumc/e. Evidently, the bulk phonons can not account 
for the resonance, since their frequency depends on q z , 
while the selection rule includes lateral momentum only. 

Our primary samples are lithographically defined Hall 
bars cleaved from modulation-doped GaAs-Alo.3Gao.7As 
heterostructures of high-mobility fi s» 3 x 10 6 cm 2 /Vs. 
The wafers are grown by molecular-beam epitaxy on the 
(001) GaAs substrate. At low T, the density of the 
2DEG, n e (in units of 10 11 cm -2 throughout the text), 
can be tuned by a combination of illumination from light- 
emitting diode and the NiCr front gate potential. The 
experiments were performed in a variable-temperature 
4 He cryostat equipped with a superconducting magnet, 
employing a standard low-frequency lock-in technique for 
resistance measurement. 

In Fig. [I] we show the normalized low- field magne- 
toresistivity p X x(B) / p xx (0) measured at T = 4 K, for 
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FIG. 1: pxx(B) I p xx (0) traces (shifted vertically for clarity) 
are shown for three densities n e of 2.05, 2.27 and 2.55 x 10 11 
cm -2 , respectively; arrows indicate the maxima for I = 
1, 2, 3, 4 and the shift of the primary (I = 1) peak with in- 
creasing n e ; Inset shows that the oscillations are periodic in 
1/B. 

the electron density n e — 2.05, 2.27, and 2.55, respec- 
tively. In addition to the damped SdH commonly seen 
in a 2DEG at this T, the traces reveal new oscillations 
that appear only at B < 0.3 T. The amplitude of the 
oscillations is about 2-3 % in these traces. 

Three aspects of the observation should be highlighted. 
First, the oscillations are roughly periodic in inverse mag- 
netic field, 1/B. The arrows next to the traces indicate 
the pxx{B{) maxima (indexed as I = 1, 2, 3, 4) in this 
oscillatory structure. In the inset we plot the order of the 
oscillations, I (and —d 2 p xx /dB 2 ), vs. 1/B for n e = 2.55 
and observe a linear dependence. Such periodic oscilla- 
tions have been seen for all n e (from ~ 1.5 to 3) studied. 
Second, with increasing n e the features shift orderly to- 
wards higher B. Finally, the oscillatory structure is ac- 
companied by a negative MR background, apparently in 
the same B range where the oscillations take place. 

We have now measured over a dozen specimen (from 
five wafers), of both the Hall bar (width from 10 /zm to 
500 fim) and the square (5 mm x 5 mm) geometries, and 
consistently observed similar oscillatory structures. On 
the other hand, the significance of the ubiquitous MR 
background remains unclear. Either negative or positive 
MR has been observed, and its strength (and even the 
sign) is largely specimen and cooling-cycle dependent. 
In the following we shall focus on the analysis of the 
oscillatory structure, in particular, its dependence on n e 
and T. 

To further quantify our results, we have performed fast 
Fourier transform (FFT) on the resistance data. As an 
example, Fig. |^ shows the FFT power spectra obtained 
from the three traces in Fig. [j] |Q. Surprisingly, such 
analysis has uncovered two frequencies, marked by A and 
B. Peak A corresponds to the main period, conforming 
to the simple fit in Fig. |l|. Peak B is somewhat weaker, 
and occurs at fs ~ 1.5/a. The shift of the doublet with 
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FIG. 2: Representative Fourier power spectra of the p xx {B) 
are shown for n e = 2.05, 2.27 and 2.55 x 10 11 cm -2 , respec- 
tively, where the doublet (see text) is marked by A, B; arrows 
indicate the shift of the slow mode with increasing n e . Inset 
shows a linear dependence of f\ (filled circles) and f% (open 
circles) vs. n e ; peak A (peak B) positions are fitted with 
sound velocity ua = 2.9 km/s (uj = 4.4 km/s). 



increasing n e is marked by three arrows for the main 
peak. The FFT data have revealed a striking linear rela- 
tion between the frequencies of oscillations and the elec- 
tron Fermi wave vector. We plot (see the inset) / 2 of 
the FFT peaks against the electron density, n e , which 
has been varied from 1.47 to 2.95 in the same specimen. 
Since fcp = y/2irn e , the observed linearity indicates that 
/ oc fcf. Such a linear dependence distinguishes the new 
oscillations from SdH, as fsdH °c k%, and is exactly what 
one expects from the phonon resonance scenario proposed 
here. As such, the oscillatory structure must be viewed 
as resulting from the resonance of the 2DEG with two 
branches of the interface modes. Using Eq. [l] and a sin- 
gle known material parameter, the GaAs band electron 
mass m w 0.068m e , we fit the data (solid line in the inset) 
and deduce a velocity for the slow (fast) mode u A « 2.9 
km/s (ub ~ 4.4 km/s). Within the experimental error of 
10% the data from several specimen collapse on the same 
lines, indicating that the new oscillations are generic in 
high-mobility 2DEG in GaAs-AUGai-^As heterostruc- 
tures. 

The T-dependence of the oscillations is consistent with 
a thermally excited phonon-scattering model. Fig. || 
shows the p xx {B) at selected temperatures (1.9 K < T < 
9.1 K ) where the evolution of the oscillations is clearly 
seen. Notice first (see inset) that p xx (0) grows linearly 
with T, indicating that acoustic-phonon scattering dom- 
inates the electron mobility in this temperature range 
H |{|. Considering the interface phonon modes of inter- 
est here, we use the value of the slow mode ua = 2.9 
km/s to estimate a characteristic temperature, T c , from 
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analog of the formula, first derived by Titeica p 



FIG. 3: p xx (B) (for n e = 2.27 x 10 11 cm -2 ) at selected T 
showing that oscillations are best developed at T = 3 — 7 
K and are strongly suppressed at higher and lower T. Inset 
shows a linear dependence of 1/fi vs. T. 



ksT c — huA^kp). The value of T c rs 5 K can quali- 
tatively account for the temperature dependence of the 
main features of the oscillations. While the SdH grad- 
ually diminishes as T increases, the oscillations are best 
developed at T m 3 — 7 K and are strongly damped at 
both higher and lower T. At T <C T c the number of 
interface phonons carrying q — 2kp becomes small and 
therefore the amplitudes diminish. At high T the smear- 
ing of the LLs prevails and the oscillations disappear as 
well. 

We now turn to the details of the theoretical explana- 
tion of the novel oscillations. We have performed calcula- 
tions ]l0| of LIP modes for the GaAs-Alo.3Gao.7As inter- 
face on the basal (001) plane. In the anisotropic case the 
speed of LIPs depends on angle between q and the [100] 
direction. Using the elastic moduli of the bulk lattices 
jnj we found a series of modes with weak anisotropy and 
a small imaginary part of the velocity (ui/ur < 0.03). 
We have studied the modes within the interval of ve- 
locities 2.4—6.0 km/s. Two close groups of modes have 
been found, one within the interval of 3—3.5 km/s and 
the other within 4.2—4.5 km/s. These modes may be re- 
sponsible for the two periods of oscillations which have 
been observed. The frequencies of the other modes found 
are too high to be detected in our experiment. Note that 
different modes may interact with electrons with different 
strengths. 

To calculate the transverse conductivity due to the 
scattering of the 2DEG by the LIPs, we employ a 2D 
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n,n' k y ,k' y 9x,<J y 

xNif n (l- f n ,)6 ky - k , y+qy 5(uj c (n' - n) - qu). (2) 

Here A is the area, Ni = [eicpituo / UbT) — 1) — , /„ = 
(exp((i5„ — fij/ksT) + i)~ , A = y/hc/eB is the mag- 
netic length, and |C(<z)| 2 = v(q)/A is the square modu- 
lus of the 2DEG-LIP interaction, which has a power law 
dependence on q. This formula can be interpreted in the 
following way. A 2D electron in a magnetic field has a 
wave function which is a product of a plane wave in the 
y direction and an oscillatory wave function, centered at 
the position xq — —chk y /eB: \& — exp(ik y y)(/) n (x — xq), 
where n is the LL index. In the absence of scattering 
the electric current may flow only in the y-direction, pro- 
viding the Hall effect. A transverse conductivity appears 
because an electron transfers wave vector q y = k' y — k y 
to a scatterer. This is equivalent to a jump in the in- 
direction at a distance Axo — chq y /eB. In Eq. this 
physics is applied to electron-interface phonon scattering. 
The mechanism of the 2DEG-LIP interaction, which may 
be either deformation potential or piezoelectric interac- 
tion, is not particularly important for our purpose. The 
overall scattering is of the same order as the bulk phonon 
scattering, since the energy densities of both excitations 
are of the same order in the vicinity of the interface. 

If the interface phonon has no attenuation, the square 
of matrix element I nn > is given by fl3|| 
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where b = qX and L n (x) is the generalized Laguerre poly- 
nomial. Substituting summation over wave vectors by 
integration in Eq. (g) one obtains 



NiUil- f n+l ) 
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dqq 3 v(q)\I nn+ i(q\)\ 5 (o<J c l - qu) . (4) 



Taking into account the imaginary part of the LIP fre- 
quency, lo = q(un — iui), we can substitute for the 5- 
function in Eq. (^) a Gaussian distribution with appro- 
priate dispersion a = qui. Since the dispersion is small 
we can set q — uj c I ju r everywhere except for the strongly 
oscillating function |/ n ;(gA)| 2 . Then after averaging we 
obtain for the transverse conductivity 



e 2 /2irh 
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FIG. 4: F 17 ,i{B) in the limit a = (dashed line) and for 
LIP (solid line). Inset illustrates the origin of the strongest 
peak in F„i(B). 

where the function F n i can be expressed as a series of 
Hcrmite polynomials of imaginary argument: 

_ {lo c I/u) 3 ( a htuj 2 \ 
~ VTT^ eXP \ l + a2mu 2 ) 

A n\(-l) 1 ( n + l\ fn + l\ 
{n + l)\k\j\ \n-k) \n-jj 



H 2 (k+j+i) \ Uy 5^7fe J 

X [2vTTq| 2 ( fc + ir H) ' 

with a = (u/a\) 2 . Hereafter, we assume that u is the 
real part of the LIP velocity. In Fig. |^ we plot F n i 
for n = 17 and I = 1 as a function of B for LIP with 
a = ll> c uj/uji (solid line) and in the limit a = (dashed 
line). As we can see, once attenuation is introduced, only 
one strong peak remains that corresponds to Eq. (|l|) at 
1 = 1. This means that only phonons with wave vector 
q = 2^/2n/X effectively interact with electrons under the 
condition n» 1. In fact, due to the Fermi distribution 
in Eq. (Q) only the values of n rs Ep/huj c are important. 
Then, indeed, 2\[2~njX = 2y/2nn e = 2kp, and we arrive 
at Eq. (|l|) for 1 = 1. The same conclusion holds for any 
I < n. 

This is an important result of our work. It can be in- 
terpreted from the following semi-classical consideration 
(see the inset in Fig. ||). Let us consider n ~ n' ^> 1. 
Since the square of the matrix element in Eq. (^|) depends 
on q only, we can put q x = 0. Then the integrand in (|^) 
is an overlap of two oscillatory wave functions shifted 
with respect to each other. In the vicinity of the turn- 
ing point the wave function always has a maximum since 
the momentum is small and the particle spends most of 
its time there. There are three possibilities. Cases 1 
and 3 in the inset show situations when turning points 
are apart from each other, and case 2 occurs when the 
turning points coincide in space. Obviously, in case 2 
the overlap integral has a maximum. This occurs when 



mu> 2 (Axo) 2 /8 = nhuj c , which is equivalent to the above 
condition qX « 2y2n. Note that the other maxima in 
Fig. U can be smeared very easily because their widths 
are proportional to n -1 / 2 , while the first maximum near 
the turning point can be approximated by an Airy func- 
tion and its width is independent of n. Thus, the maxima 
in F n i(B) for different I give rise to oscillations in p xx (B). 

As a whole, the results provide good agreement with 
our experimental data. In particular, the slow mode ua 
can be identified with the lower bunch of modes calcu- 
lated here. Within the experimental uncertainty we are 
unable to find any anisotropy for the velocity, therefore, 
the data must be viewed as an average over all directions. 
Likewise, the velocities of fast modes coincide with ub — 
4.4 km/s, but this should be taken with caution. Since 
our experiments have so far been centered on a temper- 
ature range around 5 K, a positive identification of the 
fast mode awaits for a more detailed T-dependence study 
at higher temperatures. 

In conclusion, we have discovered a new class of magne- 
to-oscillations in a high-mobility 2DEG and interpreted 
it as a magnetophonon resonance with leaky interface- 
acoustic phonons. Owning to their 2D characteristics, 
the leaky interface modes play a unique role in the scat- 
tering of 2D electrons in GaAs-Al^Gai-ajAs heterostruc- 
tures and quantum wells. This role has never been stud- 
ied before. 
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